HEAT KERNEL EXPANSIONS ON THE INTEGERS AND 
THE TODA LATTICE HIERARCHY 



PLAMEN ILIEV 

Abstract. We consider the heat equation ut = Lu where L is a second- 
order difference operator in a discrete variable n. The fundamental 
solution has an expansion in terms of the Bessel functions of imaginary 
argument. The coefficients a^(n,m) in this expansion are analogs of 
Hadamard's coefficients for the (continuous) Schrodinger operator. 

We derive an explicit formula for oik i n terms of the wave and the 
adjoint wave functions of the Toda lattice hierarchy. As a first appli- 
cation of this result, we prove that the values of these coefficients on 
the diagonals n = m and n — m + 1 define a hierarchy of differential- 
difference equations which is equivalent to the Toda lattice hierarchy. 
Using this fact and the correspondence between commutative rings of 
difference operators and algebraic curves we show that the fundamental 
solution can be summed up, giving a finite formula involving only two 
Bessel functions with polynomial coefficients in the time variable t, if 
and only if the operator L belongs to the family of bispectral operators 
constructed in f 8. . 
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1. Introduction 
The fundamental solution of the heat equation 

ut = Uxx + V(x)u, u(x, 0) = 5 y {x) 
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has an asymptotic expansion of the form 

(x-y) 2 / oo 




as t — > 0+, 



where the Hadamard's coefficients H n (x,y) are defined in a neighborhood 
of the diagonal x = y. Restricted on the diagonal, the coefficients H n (x,x) 
become differential polynomials of the potential V(x) and can be used to 
define the Korteweg-de Vries (in short KdV) hierarchy. There are several 
proofs and applications of this fact, see for instance [H [251 EH EH] and the 
references therein. One of the important steps in these papers is to use the 
connection between the heat kernel and the resolvent of the corresponding 
Schrodinger operator and the work of Gelfand and Dickey [12]. In particular, 
this construction was used to obtain explicit formulas for the higher KdV 
equations exploring different properties of Hadamard's coefficients. 

On the other hand, after the works [9l I28|. [30] . we have a much better un- 
derstanding of the KdV hierarchy and the parametrization of its solutions. 
Therefore it seems natural to reexamine the connection between the heat 
kernel and the KdV hierarchy within the framework of Sato's theory and to 
use the soliton technology as a tool to investigate the heat expansions. This 
idea was developed in [21] and it led to a new formula for the Hadamard's 
coefficients in terms of the r-function (or equivalently the wave and the ad- 
joint wave functions) of the KdV hierarchy. This formula made transparent 
some of the basic properties of the Hadamard's coefficients such as the sym- 
metry about the diagonal, or the connection between H n (x, x) and the KdV 
flows. We used this formula also in [22] to prove that heat kernel is finite 
(i.e. H n (x,y) = for n large enough) if and only if the potential V(x) is 
among the rational solutions of the KdV hierarchy (for specific values of the 
time variables) studied in [US]. These operators appear also as solutions of 
the bispectral problem |11[ 113] . For multivariable versions of the heat ker- 
nel and interesting connections with the bispectral problem and Huygens' 
principle see El [8] . 

It is natural then to ask if one can use the soliton technology to study 
heat kernel expansions where very little or even nothing is known. The 
aim of the present paper is to show that there is a very close connection 
between the heat kernel expansion on the integers and the Toda lattice 
hierarchy. Some partial results in this direction were already established in 
|14 [ 1151 [T6] , mainly in connection with the bispectral problem. Jointly with 
F. A. Griinbaum [15] we studied the heat kernel for specific second-order 
difference operators which were introduced in [18] as solutions of a difference- 
differential version of the bispectral problem. These second-order difference 
operators depend on free parameters and, appropriately normalized, they 
provide rational solutions of the Toda lattice hierarchy. The main result in 
[15] is that the fundamental solution corresponding to these operators can be 
written as a finite sum of Bessel functions of imaginary argument, which is a 
discrete analog of the finiteness property mentioned above for the continuous 
Schrodinger operator. The case when the operator L corresponds to a soliton 
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solution of the Toda lattice was considered in [16]. Both papers rely on the 
fact that the second-order operators belong to rank one commutative rings 
with spectral curves having specific singularities, which prevents their use 
for arbitrary L. The goal of this work is to derive a general formula for the 
analogs of Hadamard's coefficients in terms of the wave and the adjoint wave 
functions for the Toda lattice hierarchy and to present a few applications. 
In order to state the main results, let us first introduce some basic notations 
and define the analogs of Hadamard's coefficients in the discrete case. 

We denote by E the customary shift operator acting on functions /(n) = 
f n of a discrete variable n G Z by 

Ef(n) = f(n+l). 

For a second-order difference operator L of the form 

L = E + b n U + a n E-\ (1.1) 

the fundamental solution (or the discrete heat kernel) is the solution u(n, m; t) 
of the heat equation 

with initial condition 

u(n,m;0) = 8 n>m . (1.3) 

The papers [HI [151 05] suggest to look for a solution of (|1.2p - ([1.3|) of the 
form 

oo 

u(n, m; t) = } j a k (n, m)I n - m -k(2t), (1.4) 
fc=o 

where Ij(2t) is the Bessel function of imaginary argument. The coefficients 
ak(n,m) are the analogs of the Hadamard's coefficients H k (x,y). Plugging 
(|1.4p in (TO]) , using the identity 

d t I k {2t) = / fe+ i(2t)+/fc-i(2t) 

and comparing the coefficients of I n _ m _fc + i(2i) we get 

a k (n, m) + a fc _ 2 (n, m) = a k (n + 1, m) + b n a k -i{n, m) + a n a k _ 2 (n - 1, m), 

(1.5) 

with the convention that «j(n, m) = if i < 0. From (jl.4p and using that 
ij(0) = 5jfi we obtain 



u(n, m; 0) 




if n < m 
if n > m, 

which combined with (|1.3|) leads to 

a n _ m (n, m) = 8 n>m for every n>m. (1.6) 

It is clear that the system of equations (jl.5p - (|1.6B has a unique solution. 
Indeed, for k = equation (|1.5p gives ao( n > m ) = ao( n + ^, m )i which 
essentially means that ao(n, m) = [5 m is independent of n. On the other 
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hand, from equation (|1.6p we see that ao(rn,m) = 1, which shows that 
ao(n, m) = 1. Similarly, for k = 1 we obtain the system 



which uniquely determines ai(n, m). Next we can compute a2( n , m ), ota(n, m), 
etc. We want to stress at this point that although the system (jl.5p - (|1.6p 
has a unique solution, the coefficients in the expansion (jl.4p are uniquely 
determined from u(n, m; t) only if n < m. Indeed, for r E Nq we have 



If n < m then ak(n,m) are uniquely determined by comparing the coeffi- 
cients of the different powers of t in (|1.4|) . However, for n > m the expansion 
will contain both I r (2t) and I- r (2t) for r = 1, . . . , n — m and therefore only 
the sum of these 2 coefficients is uniquely determined by u(n,m;t). If we 
want to make the expansion unique, we can rearrange it so that the sum 
(|1.4p contains only terms Ij with j < 0. However, this would complicate sig- 
nificantly the recurrence relation (II. 5p by introducing different cases when 
n > m and n < m. Therefore we will keep the form (jl.4p with ak(n,m) 
uniquely determined by (ll.5p - (jl.6p . 

Let me now briefly describe the main results in the paper. In the next 
section we collect some basic facts about the Toda lattice hierarchy. In 
Section [3] we derive simple formulas for the coefficients a^(n, m) in terms 
of the wave and the adjoint wave function for the Toda lattice hierarchy 
(or equivalently the r- function). In other words, the soliton technology 
allows us to integrate explicitly equations (jl.5p - (ll.6p and write ak(n,m) in 
closed form. As a first application of this formula, we show in Section 0] 
that «fc(n, n) and ctk+i(n,n — 1) generate a hierarchy which is equivalent 
to the Toda lattice hierarchy (notice that we need two coefficients in order 
to generate two equations for a n and b n ). This provides a discrete analog 
of the connection between the Hadamard's coefficients H n (x,x) and the 
KdV hierarchy. It would be interesting to see if some of the proofs in the 
continuous case can be adapted to this situation. We note, however, that 
the situation in the discrete case is a little bit more complicated because 
the fe-th vector field of the Toda lattice X/j is not equal to the vector field 
X' fc generated by ctfc(n, n) and ak+i(n, n — 1), but we have X' fc = X^. + 
a linear combination of ~Kk-2j with 2j < k, see Theorem l4.4i All this follows 
easily from the explicit formula for afc(n,m), but clearly a direct approach 
should be carefully constructed in order to capture the lower-order vector 
fields that appear. 

Next we prove that the heat kernel can be written as 



ai(n, in) = a\(n + 1, m) + 6 ; 
ai(m + 1, m) = 0, 



71 




(1.7) 



(1 +pi(n, m;t))I n - m (2t) +p 2 {n,m]t)I 1 



i(2t) 



(1.8) 



where pi(n,m;t) and p2(n,m;t) are polynomials in t with coefficients de- 
pending on n and m, such that p\{n, m; 0) = p2(n,m;0) = if and only 
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if the operator L belongs to the family of bispectral operators constructed 
in [18] . The "if" part is essentially the main result in [15] (see Section [6] 
for more details). In order to prove the "only if" part, we first discuss in 
Section [5] the operators constructed in [18] and we show that they can be 
characterized by the vanishing of a specific linear combination of the Toda 
flows, after a particular point. The heart of the proof is to show that if 
this specific linear combination of the Toda flows vanishes after a particular 
point, then the operator L belongs to a rank one commutative ring with 
spectral curve of the form v 2 = (u — 2) 2Nl+1 (u + 2) 2N2+1 , which leads to the 
operators introduced in [18]. Using this result and the connection between 
the heat coefficients and the Toda flows, we establish the "only if" part in 
Subsection 16.21 

The results in the last section deepen the mystery surrounding the bis- 
pectral problem and its connection with the heat equation. The finiteness 
property of the heat kernel provides one more reason to believe that the 
operators constructed in [18] are the discrete analogs of the Adler-Moser op- 
erators pp. It would be interesting to see if some of the recent developments 
|19[ [20] in the purely continuous version of the bispectral problem can be 
extended to the operators in [18], or to the more general (higher-order) bis- 
pectral operators in [T7], that we constructed in collaboration with L. Haine, 
inspired by Wilson's work |34[ I35j. 

2. The Toda lattice hierarchy 

In this section we briefly recall the construction of the Toda lattice hier- 
archy, its wave and adjoint wave functions. For more details we refer the 
reader to the paper [31] , or to the more recent contributions [21 [32] . 

Let us first introduce some notations, which will be used throughout the 
paper. We denote by A and V the customary forward and backward differ- 
ence operators, acting on a function f(n) = f n as follows 

A/(n) = f(n + 1) - f(n) = {E — Id)/(n) 

V/(n) = f(n) - f(n - 1) = (Id - E' 1 )/^). 

Sometimes, we will write A n and V n if the operators are applied to functions 
depending on several variables to indicate that A and V act in the variable 
n. 

When we work with pseudo-difference operators of the form 
M = EL-oo c k (n)E k we will denote by M+ = Y%=o c k{n)E k (resp. M_ = 
X^fe=-oo c k(n)E k ) the positive (resp. negative) difference part of M. 

For a second-order difference operator L = E + b n ld + a n E~ l , the Toda 
lattice hierarchy is defined by the Lax equations 

Fit 

— = [(Li)+,L], for j = 1,2,.... (2.1) 

It is well-known that the vector fields Xj(L) = [(L J ) + , L] commute with each 
other, that is, each of the equations (|2. If) defines a symmetry of any other 
equation. For this reason, the family of equations (|2.1|) is called a hierarchy. 
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The first equation, corresponding to j = 1, is the well-known Toda lattice 
equation. 

We denote by W n = 1 + YlT=i' l l ; k(n)E~ k the wave operator, i.e. the 
operator which conjugates L to E 

L = W n EW-\ (2.2) 

The flows (|2.ip can be extended on W n by 

dW n 



dSj 



-{U>)-W n . (2.3) 



The wave function ^ n (z;s) and the adjoint wave function ^*(z;s) are 
defined by the formulas 



oo 

■j 



\i =1 

OO , / n \ / oo 



and 



z 

k=l / \j=l 



oo 

■3 



K(z; S ) = (w-l 1 yz- n e W [-Y J t 



i=i 



where = E~ 1 . With these definitions we have 

L*„(z; = s) and ^ ' ; = (L 3 )+\I> w (z; 8 ). (2.5) 

We denote by ^f n (z; s) and ^* n {z; s) the reduced wave function and the 
reduced adjoint wave function obtained from ^ n (z;s) and ^*(z;s), respec- 
tively, by omitting the exponential factors, i.e. 

^ ;s) = 1 + £Mj£> and r .(. ; .)-i + ±^p.. (2.6) 

k=l k=l 

In order to simplify the notations, we will often omit the s-dependence, 
e.g. we will write simply ^> n (z) instead of ^ n (z;s), etc. Notice that the 
first equation in (|2.5p is equivalent to the following equation for the reduced 
wave function 

zV n+1 (z) + bjf n (z) + a n ^ n ' 1 ^ = zV n (z). (2.7) 



oo 
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For series CkZ k and for formal pseudo-difference operators dkE k 
we denote 




res z y CkZ = C-i and resg 



It is easy to see that for pseudo-difference operators P n = ^ fc dk{n)E and 
Qn = Ylk c k(n)E k we have 

res, ((P n z n )(Q* n _ lZ ' n )) = ms E {P n Q n ). (2.8) 

From this equality it follows that ^ n (z; s) and ^* n {z; s) satisfy the following 
bilinear identities 

xea^ffidfr ■■■8$ i * r + ko (z i a))** n (x-,aj) = 0, (2.9) 



where ko, k\, . . . , kj G No and dj = Indeed, from the second equation 
in (|2.5|) we see that it is enough to prove (|2.9|) for k\ = k<z = ■ ■ ■ = kj = 0. 
Using now ()2.8[) we get 



res z (* n+ko (z;s))*Uz;s)) = res, ([E^W n z n ){{W~\)* z~ n )) 
= ies E {E ko ) = 0, 



which establishes (j2.9j) . 

For example, for /co = and ko = 1 with /cj = for j > 1, the identity 
(1231) gives 



^i(n) + ^f*(n) = (2.10a) 
V> 2 (n + 1) + Vi(n + 1)^1 (n) + $j(n) = 0. (2.10b) 

From (12. 2h and (|2.4() we can express 6 n and a n in terms of the coefficients 
of ^ n (z; s) and ^* n (z; s) 

b n = ip 1 (n) + ipt(n + l) (2.11a) 
«n = Mn) + ^i(n)^*(n) + ip 2 (n). (2.11b) 

Remark 2.1. If we fix n and take k$ = 0, then the bilinear identities 
(|2.9p show that ^> n (z; s)z~ n and ^(z; s)£ n are the wave and the adjoint 
wave functions of the Kadomtsev-Petviashvili (KP) hierarchy. Thus, by the 
classical theory (see [9], [TOj Theorem 6.3.8, p. 97]), there exits a r-function 
T n {s) such that 

*n(*;.) = T " (S " [ f 1]) and K(z;s) = ^ 

where [z] = (z, z 2 /2, z 3 /3, . . . ). This function plays a central role in the KP 
and the Toda lattice hierarchies. Since we are not going to make an explicit 
use of it in the present paper, we stop the discussion here. Notice that the 
last formulas allow us to express the coefficients of *ff n (z; s) and ^(z; s) in 
terms of r n (s). In particular, this means that all formulas involving ^> n (z; s) 
and ^n(z; s) that follow can be rewritten as formulas involving only T n (s). 
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3. Explicit formulas for the heat coefficients ak(n,m) 

In this section we show that, using the notations in the previous section, 
we can "integrate" equation (|1.5j) for all k £ N and obtain simple formulas 
for ak(n,m) in terms of the wave and the adjoint wave functions of the 
Toda lattice hierarchy. Before we present and prove the general formula, we 
illustrate how the method works for ot\{n,m) and a 2 (n,m). This will help 
the reader understand the nature of these formulas and the importance of 
the bilinear identity (|2.9j) . 

Plugging k = 1 in (jl.5p and using (|2.1ip we get 

«i(n, m) = «i(n + 1, m) + ipi(n) + ipl(n + 1). 

Using (|2.10ap we can rewrite the last equation as 

ai(n, m) — V'i(ra) = «i(n + 1, m) — ip\(n + 1), 

which shows that a\(n, m) — tpi(n) = (3i{m) is a function independent of n. 
From (jl.6p we see that ai(m + l,m) = 0, which gives f3\(m) = —ip{m + \) = 
ip*(m + 1). Thus we have 

ai(n,m) = ^i(n) + ^{{m + 1) = res* [^ n (z)^* m+1 {z)] . (3.1) 

Similarly, for k = 2 we obtain the following equation for a 2 (n, m) 

a 2 (n,m) + 1 = a 2 (n + l,m) + [ipi(n) + tp$(n + l)][V>i(n) + V4( n + 1)] 

+ Mn) + Mn)^i(n) + ^{n). (3.2) 

Using now both equations in (|2.10p one can check that (13. 2p is equivalent to 
the equation 

A n (a 2 (n, m) - fain) - ^i(n)^*(m + 1) — n) = 0, 

which means that a 2 (n,m) = ^(n) + ^i(n)^ 1 (m + 1) + n + /3 2 (w.), where 
/3 2 (m) depends only on m. From (|1.6p we see that a 2 (m + 2, m) = leading 
to /3 2 (m) = -^ 2 ("i + 2) —ipi(m + 2)'4>l(m+ 1) - m-2 = -0 2 ("i + 1) - m - 2, 
where in the last equality we used again (|2.10bp . Thus we obtain 

a 2 (n, m) = fa(n) + il)i{n)ij)\{m + 1) + ip^i 171 + 1) + (n — m — 2) 



= res, p + n-rn-2) j^±lM j. (3 . 3 ) 

The point is that this process can be successfully iterated providing simple 
formulas for ak(n,m) for every k £ N. Below we introduce the necessary 
functions which are needed to state and prove the general formula, extending 
(13TB and D. 



Let us define monic polynomials Q^.(z) by the formula 



fk 
fc-i 



fl2M-«*+0-*)E (7-7- 1) ft (3 - 4) 

From (|3.4p it is easy to see that 

<rfw-««e« = nr*(£:l) < 3 - 5 > 
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and 



OfW + QU*) - *Qt\{z) = (/ ■ (3-6) 



Let us denote 



Q n k J™(z 2 ) if k is even 
«Q(P5^(* 2 ) if k is odd. 



< 7fc (n,m;z) = <(_^; i ;_ 2x """^T (3-7) 



It is clear that g k (n,m; z) is a monic polynomial in z of degree k, and it is 
an even polynomial when k is even and an odd polynomial when k is odd. 
From the defining relation (|3.7p and equation (|3.5p one sees immediately 
that 



gk(n,m; z) — zgk~i(n — l,m; z) = if A; is odd (3.8a) 



and 



g k (n,m;z) - zgk-i(n - l,m; z) = g k (n,m; z) - z 2 g k - 2 {n - 2,m;z) 

2(n — m — k) (n — m — 1\ „ , . 

1 ] if k is even. (3.8b) 



k \ k/2-1 

Now we are ready to formulate the main result in this section which gives 
an explicit formula for a&(ra, m) in terms of ^? n (z) and ^*(z). 

Theorem 3.1. The coefficients ak(n,m) in the expansion (jl.4[) of the fun- 
damental solution of the discrete heat equation (|1.2|) can 6e expressed in 
terms of the reduced wave function ^ n {z) and the reduced adjoint wave func- 
tion ^* n {z) as follows 



a k {n,m 



res. 



g k (n,m;z) 



TO+l 



2 



(3.9) 



where g k (n,m; z) is the polynomial defined by (|3.4p and (|3.7p . 

Proof. We need to check that the functions defined by formula (|3.9p satisfy 
the difference equation ()1.5p and the initial conditions (|1.6p . Let us first 
establish that equation (]1.5p holds. We denote by RHS (resp. LHS) the 
right-hand side (resp. the left-hand side) of equation (|1.5p . Plugging ()3.9p 
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in (|1.5j) leads to 

RHS = a k (n + 1, m) + b n ak-i(n, m) + a n a k - 2 (n - 1, m) 



res. 



g k (n + l,m;z) h b n g k _ 1 {n,m; z) 



+ a n 9k-2{n- l,m;z) 



^n-l{zW m+1 {z 



res. 



((z 2 - zb n )g k - 2 {n - 1, m; z) + b n g k -i(n, m; z)) 



tf„(z)^ +1 (z) 



+ (flfcC" + 1, m; z) - z 2 g k _ 2 (n - 1, m; z)) 



tf n+ i(z)#* x (z) 



where in the last equality we eliminated \P n _i(z) using (|2.7p . 

Now we consider two cases depending on whether A; is even or odd. 

Case 1. Let k be even (hence A; — 1 is odd). Using f)3.8f) we see that 

(z 2 - zb n )g k ^ 2 (n - 1, m; z) + b n g k ^i{n, m; z) = z 2 g k - 2 (n - 1, m; z) 

and 



+ 1, m; z) - z 2 g k - 2 {n - 1, m; z) 



2(n + 1 — m — A;) / n — m 



k 



[k/2 



1 



which lead to the following expressions for RHS and LHS: 

$ n (z)V* m+1 (z) 



RHS = res, 



z g k - 2 (n - l,m; z)- 



2(n + 1 — m — k) / n — m 
+ k \k/2 - 1 



and 



LHS = res. 



(g k (n, m; z) + g k - 2 (n, m; z)) 



Using now equation f)3.6|> we get 



g k (n, m; z) + g k - 2 (n, m; z) - z 2 g k - 2 (n - 1, m; z) 
= Q n k J?{z 2 ) + Qlj^iz 2 ) - z 2 Qlj™-\z 2 ) 
2(n + 1 — m — k) ( n — m 



k 



\k/2-l)' 



which shows that RHS = LHS and establishes (|1 .5f) . 
Case 2. Let k be odd (hence A; — 1 is even). 
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Using again (13. 8p we obtain 

g k (n + 1, m; z) - z 2 g k _ 2 {n - 1, m; z) 

= z(9k-l(n, m; z) - z 2 g k _ 3 (n - 2, m; z)) 
2(n — m — fc + 1) / n — m — 1 
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fc-1 



(fc-l)/2-l 



and 



(* 2 



zb n )g k - 2 (n - l,m;z) + b n g k -i(n, m; z) 
= (z 2 - zb n )zg k ^ 3 {n - 2, m; z) + 6 n # fc _i(n, m; z) 
= b n (g k -i{n,m; z) - z 2 g k - 3 (n - 2,m;z)) + z 3 g k - 3 (n - 2,m;z) 



2(n — m — fc + 1) / n — m — 1 



Thus 

RHS = res. 



V0- 1 )/ 2 - 1 

2(n — m — fc + 1) / n — m — 1 



+ z 3 g k - 3 (n - 2,m;z). 



1 



+ z 3 g k - 3 (n - 2,m;z) 



, (*¥„+i (z) + 6 n ^ n (z))^±i^ 

(fc - l)/2 - + V ; -r n nv . 



res. 



2(n — m — fc + 1) / n — m — 1 
V(fc - l)/2 - 1 

X ^n(z)^„ +1 (z 



+ z 2 g k _ 3 (n - 2,m;z) 



where in the last equality we used (|2.7|) to eliminate z^ n+ \{z) + b n ^ n (z). 
On the other hand we have 

LHS = res^[(5 fc _i(n - l,m;z) + g k - 3 {n - 1, m; z))^ n (z)^ +1 (z)] , 

and applying the definition (|3.7[) and equation (|3.6j) we get 

g k -i{n - 1, m; z) + #fc_ 3 (rt - 1, m; z) - z 2 g k _ 3 (n - 2, m; z) 

— ^(A,-1)/2V Z / + V(H)/2-ll z J z V( fe _i)/2-ll z J 

2(n — m — fc + 1) / n — m — 1 \ 

= fc^l V(fc-l)/2- lj' 

which shows that RHS = LHS and completes the proof of the difference 
equation (jl.5p . 

It remains to check that the initial condition (jl.6l) holds, i.e. we need to 
show that for every k £ N, a k (m + fc, fc) = 0. Using the defining relations 
(|3.4p and (|3.7p one can easily see that 

g k (m + fc, m; z) = z fc . 

Thus (pT9]) gives 



a k {m + k,m) = res 2 



res z [* m+fc (z)^ +1 (z)] =0, 



where in the last equality we used the bilinear identity 



□ 
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4. Generating the Toda flows with the heat coefficients 



As a first application of formula (|3.9[) we prove in this section that 
A n ak+i(n, n— 1) and a n V n oek(n, n) generate the Toda lattice hierarchy. This 
is a discrete analog of the remarkable connection between the restriction 
of the Hadamard's coefficients on the diagonal and the Korteweg-de Vries 
hierarchy. Before we prove this, we establish some auxiliary facts. 



Proposition 4.1. Let us denote 

t*(n) = res 2 (jv n (z)%(z) 
l k (n) = res z (V 5 " 1 *^ (;<:)#; (z)J . 

Then the Toda lattice hierarchy (|2.1|) is equivalent to the equations 

t k (n+l) -t fc (n) 
: a n {\ k (n + 1) - [fc(n)). 



ds k 
da n 



ds h 

Proof. Notice that equation (|2.4bp implies 

oo 

w- l = i + J2 E ~ j ■iPK n + 1 ) 

and therefore, using (|2.2p and (|2.4a|) . we get 



(4.1a) 
(4.1b) 

(4.2a) 
(4.2b) 



L k = W n E k W~ x = J2 



m=0 



J, ip m -j(n)ipj(n + 1 + k 

3=0 



m) 



E 



k—m 



where ipo( n ) = V ; o( n ) = 1- F rom the above equation it is clear that the coef- 
ficients of E° and E~ 1 in the operator L k are ^(n+1) and tfc(n) respectively. 
This shows that 

[(L k ) + , L] = [L, (L fc )_] = [E + Md + a n E-\x k {n)E- 1 + 0(E- 2 )} 
= (r i; (n + l)-t fc (n))Id + 0(^- 1 ), 

which gives (I4.2ap . Similar computation shows that the coefficient of E^ 1 
in [(L k ) + ,L] is a n (lk(n + 1) — l k (n)), completing the proof. □ 

Lemma 4.2. Lei 6 No and j € Z 6e fixed. TTien 

a k (n,n+j) £ C[a n ,b n ,a n ±±,b n ±i,a n ± 2 ,b n ±2, . . .]. In other words, 

a k {n,n + j) is a polynomial of finitely many of {a n , b n , a n ±i, b n ±\, . . . }. 

Proof. The proof can be easily obtained by induction on k and + using 
(fT3D-(fL6D. □ 
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Example 4.3. We list the values of the first few heat coefficients in a 
neighborhood of the diagonal n = m 

Qi(n,n) = b n , 

a\{n, n — 1) = 0, 

ai(n,n+ 1) = b n + b n+1 

a 2 (n, n) = a n+i + a n + b 2 n - 2, 

a 2 (n,n - 1) = a n — 1, 

a 2 (n, n + f ) = a n + a n+ i + a n+2 + 6^ + b 2 n+1 + 6„6„+i - 3. 

The main result in this section is the following theorem. 

Theorem 4.4. The system of differential- difference equations 
db n 

^p- = a k+1 (n + 1, n) - «fc+i(n, n - 1) = A n a k+ i(n, n - I) (4.3a) 

— — = a n (a k (n, n) - a k {n - l,n - 1)) = a n V n a fc (n, n), (4.3b) 

where k £ N forms a hierarchy. Moreover, if denote by X& and Xl £/te vector 
fields corresponding to the flows d/ds k and d/ds' k given by (12, ip and f)4.3|) 
respectively, then 

i=0 v 7 

where [x\ denotes the greatest integer less than or equal to x. 

Remark 4.5. Notice that according to Lemma 14. 2\ for every k £ N the 
right-hand sides of f)4.3|) are polynomials of finitely many of {a n+ j, b n+ j}j e z, 
i.e. (|4.3p is a well-defined system of differential equations for a n and b n . 
The above theorem essentially says that the hierarchy of equations (|4.3p 
is equivalent to the Toda lattice hierarchy (|2.f p modulo a simple (linear) 
change of variables given by (|4.4p . Below we write also the explicit linear 
combination that gives X^ in terms of {X^}, which will be needed later. 
First we show that for every k £ N and m £ No satisfying m < \_(k — l)/2j 
the following identity 
m 



i=0 



k — 2i\ i J \m — i 



holds. Indeed, if m = then (|4.5p is obvious. For m > 1 we can rewrite the 
left-hand side of (14.51) as 



k v-^ , s i (k — 1 — i \ ( m\ k (k — 1 — m\ 

> (- 1 ) i • = — , = 0. 

z — ' \ m — 1 } \i ) m\ k — m I 



m 

i=0 



In the first equality we used the well-known binomial identity 

— i\ f Tn\ fk — m 
r J \ i J \ k — r 



i=0 
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which can be easily proved by applying the principle of inclusion and exclu- 
sion to the following problem: In how many ways one can select r of given 
k distinct objects, so that each selection includes some particular m of the 
k objects. 

Combining (14.41) and (14.5(1 one can deduce that 

L(*-1)/2J m 

X *= £ • X U- (4-7) 
i=o v/ 

Proof of Theorem \4-4\ It is enough to prove (|4.4p . because this formula and 
the fact that {X&} commute will imply that {Xl} commute, i.e. the equa- 
tions (|4.3p form a hierarchy. 

Using (|3.4p and (|3.7p one can easily check that 

L(fc+1)/2J , ^ i\ 

S^fon-l;*)^^* g ^ ( ^ i ) ( 48 ) 
which combined with (|3.9p and (|4.1ap gives 

a fc+ i(n,n - 1) = t fc (ra) + fe — : — ( _i ) r fc-2i(™)- (4.9) 

i=l * \ « / 

On the other hand, from equation (|3.8[) it follows that 

g k +i{n,n- l;z) = zg k {n - l,n - 1; z) - 6k,i, 

which shows that 

L(fc+1)/2J , f\i /i ■ <\ c 

g k (n-l,n-l;z)=z k + k £ til ( fc 7 * " 1 ) + £M, (4 .10) 

i=l x 7 

and therefore, using ([3T9"P and (HTTbT) we get 

L(fc+i)/2j ^ . ._ . 
a fe (n- l,n-l) = I fc (n) + A; £ i-^-f k- 2i (n). (4.11) 

i=i 

The proof now follows from (|4.9p , (|4.1ip , Proposition 14.11 and the fact that 
*o(n) = 0) t_i(n) = 1, (o(ra) = 1, Li(n) = are independent of n. □ 

5. DARBOUX TRANSFORMATIONS FROM Lq = E + E^ 1 AT THE END 
POINTS OF THE SPECTRUM 

In this section we focus on certain second-order difference operators Ljvi,JVa > 
which were introduced in [IB] in connections with a difference-differential 
version of the bispectral problem [11] . These operators can be defined by 
successive Darboux transformations from the operator Lq = E + E^ 1 at the 
end points ±2 of the spectrum. Recall (see [21]) that the Darboux trans- 
formation of a second-order operator L at point cq consists of factorizing 
L — cold as a product of first-order operators and producing a new operator 
L by exchanging the factors, i.e. if we write L — cqH = VQ, then L is 
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defined by L — cold = QV. The main result in this section is a character- 
ization of these operators in terms of the vector fields of the Toda lattice 
hierarchy. This is needed in the next section where we prove that the heat 
kernel expansion for these operators can be written as a sum of only two 
Bessel functions with polynomial coefficients (in the time variable t), and 
that this property completely characterizes the operators LjVi,7V 2 - 

5.1. Constructing the operators L^ lt jsf 2 . The operators £jVi,JVa are ob- 
tained by the following sequence of the Darboux transformations 

L - 2 Id = VqQo L lj0 - 2 Id = Q V = ViQi rx ■ ■ ■ 
Ln u0 - 2 Id = QjVj-iTVj-i 

Lni,o + 2 Id = Vni Qn 1 r\ L Nlj i + 2 Id = Qn 1 Vn 1 = V^+iQ^+i r\ ■ ■ ■ 
^> Ln 1 ,n 2 + 2 Id = Qn 1 +n 2 -iPn 1 +n 2 -i- (5-1) 
At each step, the factorization of the operator Lj x j 2 ± 2Id depends on one 
free parameter. Thus, the operator Ln 1: n 2 w iH depend on N\ + N 2 free 
parameters. The operator Ln 1) n 2 belongs to a rank one commutative ring 
of difference operators, i.e. we can apply the correspondence between com- 
mutative rings of difference operators and algebraic curves developed in the 
papers [231 ISSJ [33] . We first sketch the main steps of this construction with 
an emphasize on the operators obtained by the Darboux process (|5.ip and 
refer the reader to [T71 [18] for more details. 

Following [26], we call a difference operator M = Y2k=K Hk{ n )E k prop- 
erly bordered if /Uif_(n) 7^ and HK + { n ) 7^ for all n £ Z; the inter- 
val [K_,K + ] is the support of M. A commutative ring A of difference 
operators is called rank one, if it contains two properly bordered differ- 
ence operators M' and M" with supports [K'_,K' + ] and [K!L,K'!] such that 
gcd(K'_,K'l) = 1, gcd(i^,i^) = 1 and K'_K\ < K' + K'L. In that case, 
Spec(^4) is an irreducible complex affine curve that completes by adding two 
nonsingular points at infinity. 

Starting with a properly bordered second-order difference operator L = 
E + 6 n Id + a n E~ x we denote by Al the ring of all difference operators 
commuting with L, i.e. 

f K+ 1 

A L = I M = Vk{n)E k : [M,L] =0>. 

( k=K- 

One can show that Al is in fact a commutative ring consisting of properly 
bordered difference operators and it is a rank one ring if and only Al con- 
tains an operator which is not a polynomial of L. For every M € Al the 
operators L and M satisfy an algebraic relations of the form 

M 2 = MTi(L) + T 2 (L), 

where T\ and T2 are some polynomials. This equation defines an affine 
curve fL,M( u ' v ) = v 2 — vT\{u) — ?2(m) = 0. It is easy to see that for 
every M S Al there exists a unique polynomial gu such that <?a/(0) = 
and the operator M — guij^) contains only nonnegative powers of E, i.e. 
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M = gu{L) + Y^j=P c k( n )E k . From this, it follows that the ring Al is 
generated by two operators {L, M} where M is chosen so that v{M) > is 
minimal. The spectral curve Spec(^4i) is Jl,m(u,v) = and the complete 
curve is Xl = Spec(AL) U Q^o}- 

The Baker function for Al is the unique (up to a factor independent of 
n) eigenfunction for the operators from Al- If we denote by Al the ring of 
functions meromorphic on X with poles only at Q^, then for every M £ Al 
we have 

M* n (P) = h M (P)* n (P), where h M (P) £ A L . 
Moreover, if the support of M is [-fT_, -fT_|_] then Iim(P) has poles of orders 
K + and K_ at Q+ and Q^, respectively. 

Let us denote by A^,n 2 = Al n N2 the ring of all difference operators, 
commuting with LN lt N 2 - I n [18] it was shown that AN lt N 2 is a rank one 
commutative ring of difference operators, which is isomorphic to the ring 

A NltN2 = Al Ni _ N2 = C[x + x' 1 , f Nl:N . 2 (x)] C CfoaT 1 ], 

where 

tN lt N a {X)- x jV 1+ jV 2 +l • ( 5 - 2 ) 

In other words, if ^f n is the Baker function for A^ n 2 then 

L Nl)N2 ^ n {x) = (x + x-^nix) (5.3a) 

M NuN2 ^ n (x) = f Nl>N2 {x)V n (x), (5.3b) 

for some Mn 1s n 2 £ AN lt N 2 an d the ring Ani,n 2 is generated by Ln x ,n 2 and 
Mwi ,Ar 2 • The spectral curve is given by the equation 

Spec^jv^ivJ : f 2 = (u - 2) 2iVl+1 (u + 2) 2Af2+1 . (5.4) 

Clearly, Spec^jVi.A^) is rational and has a cusp at u = 2 (resp. ii = —2) 
when iVi > (resp. N2 > 0). Another property, that we will need later, is 
that for every j £ Z we have 

x j f NuN2 (x) £ A NltN2 , (5.5) 

see [181 the proof of Theorem 4.2.] (note that x here corresponds to z + 1 
in [TBI). 

Conversely, let L = .E+^Id+an-E 1 1 be a properly bordered second-order 
difference operator whose spectral curve is given by equation (|5.4p . From 
the correspondence in we can conclude that Al is obtained from An x ,n 2 
for a specific choice of the free parameters in the Darboux steps, up to a 
conjugation by a nonzero function g n , i.e. 

A L = {g^Mgn ■ M £ ^tvJ. (5.6) 

If L is nonconstant (i.e. at least one of the functions a n and b n is not a 
constant), then it is easy to see that all second-order operators operators in 
Al with support [—1,1] must have the form 71 L + 72 Id for some constants 
71)72- From this, it follows that the only possible functions g n in (|5.6p are 
g n = c n , and we must have 

L = C - n ~ 1 L NlN2 c n + dld, (5.7) 
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for some constants c 7^ and d. 

One possible way to eliminate this freedom and characterize precisely the 
operators -LjVi,Ar 2 is to consider properly bordered second-order difference 
operators of the form L = E + 6 n Id + a n E~ l with spectral curve given in 
(|5.4p and coefficients satisfying 

lim b n = 0, lim a n = 1. (5-8) 

n^oo n— >oo 

Indeed, if L = Ln ± ^ 2 then the coefficients 6 n , a n can be computed from the 
formulas 

3 r„ +1 (s) f n+ i(s)f n _i(s) 
b n = ^— log _ , ; , a n = —-2 , 

where r n (s) is a polynomial in n, which makes (|5.8p obvious. We note, how- 
ever, that f n (s) in the last formulas differs from the r-function introduced 
in Section [2] (due to the different approach in the papers [T71 HI]). Con- 
versely, if L is nonconstant coefficient operator which belongs to a rank one 
commutative ring of difference operators whose spectral curve is given by 
(|5.4p . then (|5.7p must hold which combined with (|5.8p shows that c = ±1 
and d = 0. Conjugating L^ lt N 2 by (— l) n essentially exchanges the roles of 
+2 and —2 in (|5.ip (or, equivalently, the roles of Ni and N2). Thus, we have 
L = Liv^jVj or L = Ljv 2j tvd completing the proof in this case. Finally, if a ra 
and 6 n are constants, then (|5.8p implies that L = Lq. 

Next, we establish several new facts needed for the characterization of the 
operators L/Vi,jV 2 i n terms of the Toda vector fields proved at the end of this 
section. 

Lemma 5.1. Let A be a ring of Laurent polynomials in x such that 
An 1s n 2 C A for some N\, N2 G No- If A contains a polynomial p(x) such 
that < deg(p(x)) < N% + N2 then one of the following must hold 

(i) Ni > 1 and f Nl _i iN2 (x) G A; 

(ii) N 2 > 1 and f NliN2 ^(x) G A. 

Proof. Assume that p(x) G ^4 n C[x] is such that < deg(p(x)) < N% + iV2. 
Since x + x _1 G ^4, it is clear that p{x) +p(l/x) G A, which combined with 
p{x) G A shows that f(x) = p{x) —p(l/x) G A. Notice that f(l/x) = —f{x) 
and therefore we can write f(x) in the form 

/(*)=(*-) /(*), (5.9) 

where f(x) G C[x,x _1 ] is not divisible by x — 1/x (i.e. /(— 1) / or 
/(l) ^ 0) and 2/ + 1 < iVi + iV 2 . 

Let us denote T = mm(Ni,N 2 ) and T = max(Ni,N 2 ) > 0. Choosing 
e = 1 if T' = N2 and e = — 1 otherwise we can rewrite fNi,N 2 ( x ) as 

(x + e) 2T ' +1 (x-e) 2T+1 



fN u N 2 (x) 



X T'+T+1 
1\ 2T+1 \ 1 

x + 2e + - 

X X 



T'-T 
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Below we consider separately the cases I < T — 1 and I > T. 

Case 1: I < T — 1, Let us write f(x) in (|5.9p as follows 

5 



3=0 

Case i.a. Assume first that /o / and T' > T. Then, we can multiply 

fix) by 

1 \2(T-/) , 1 ,T'-T-l 

x ) ( x + 2e + - I £ A 

XI \ XI 



and we get 



3=0 



Notice that for j > 1 the terms in the above sum are multiples of /ni,n 2 
and therefore belong to A by (|5.5|) . Thus, the term corresponding to j = 
also belongs to A and since /o / we obtain 

x-- x + 2e + - =^ * =A= ^ 

X J \ X J X 1 +1 

which is what we wanted to show. 

Case l.b. Assume now that /q / but T' = T (hence e = 1), and let us 
multiply f(x) by 

/ l\2(T-i-l) / r 

fx ] (x-2+- ] eA. 

We have 



x / \ x 



m , x 2T-1 / 1 \ ■? / 1 

X+2+- *-2+-) £ A 

3=0 



Again all terms for j > 1 in the sum above are in ^4 and therefore, since 
/o / 0, we get 

ar - - [x-2 + -)= - - ^ - G A. 

X J \ X J X Z1 

Case I.e. Let now Jq = 0, i.e. /(— e) = 0. This implies that /(e) 7^ 0, 
because otherwise (x — l/x) will divide /(x) contrary to our factorization in 
(|5.9p . Thus we can write /(x) as 

m / 1 \ j 

f(x) = ^2f j [x-2e + -) , (5.10) 

3=0 V XJ 
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with jo 7^ 0. Using now (|5.10p and multiplying f(x) by 



1\ T - W / 1\ T '-' k 

2e + - x + 2e + - G A 

x \ x I 



we get 



m ( i\ 2l + 1 / 1 \j+T-i-i / 1\ T ' _/ 

g/i(*--J (*-* + ;) (* + & + i) ^ 

For j > 1, the j'-th term in the above sum is 

1 \ J * -1 (x - e) 2T+1 (x + e ) 2T ' +1 
x ~ 2e + x) X T+T'+1 G ^' 

and therefore for j = we have 

(x-e) 2T - 1 (x + e) 2T '+ 1 



x T+T ' 



G A. 



Case 2. Finally, if I > T we can write f{x) as 

i 2T+1 5 

f(x) = ( X 

3=3o 



l\ 2T+1 s ( IV 
-) I>(* + 2 £ + -), < 5 '"> 



with h jo / and S 1 < T'-T-l. Multiplying (foTTTT) by (x + 2e + ~) 
^4 we obtain as before that 

( x _ e )2T+l (2 , + e) 2T'-l 



IxT'-T-l-jo 



X T+T ' 



G A 



which completes the proof of the lemma. □ 

Corollary 5.2. If A is a ring of Laurent polynomials in x, such that 
-Ajv.jv G ^4 /or some iV G No, f/ien ^4 = An u n 3 for some Ni,N% G No- 

Proof. For every /(x) G A there exists a unique polynomial p(x) such that 
p(0) = and f(x) — p(x) G C[x + x -1 ]. Since x + x" 1 £ A we see that A 
is generated by x + x~ l and the polynomial of minimal positive degree in 
AflC[4 The proof now follows immediately from Lemma 15,11 □ 

Let us denote by q k {x) the odd polynomial of degree 2k + 1 given by 

k 

%(^=Er/ 1 )(- i )^ 2fe_2j ' +1 - ( s - 12 ) 



3=0 

Then 



I \ 2fc+l 



X 



fk,k( x ) = [ x = q k (x) - q k (l/x). (5.13) 



If we denote u = x + x 1 then from (|5.12j) and (|5.13j) it is easy to see that 

^-A)a k (x) + (-l) k+1 ( 2k + 1 



q k+1 (x) = (u 2 -4:)q k (x) + (-l) k+L ^ fc \u. (5.14) 

Using the last relation, one can deduce by induction on k that 

q k (x) + gfc(l/x) = P k (u), (5.15) 
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where 

k{ } ^ j\ (2fc + l-2j)!! 

(2fc + l)! (fc-j)! 2fc _ 2?+1 

fc! ^ j! (2fc + l-2j)! 

and (2j + l)!! = 1 ■ 3 • • • (2j + 1). 

Proposition 5.3. Let Lw lt N 2 ^ e ^ e second-order difference operator con- 
structed in ()5,ip and A?" = max(Ai, N2). Then for every k > N we have 

Pk{LN 1 ,N 2 )+ ^ An ± ,N 2 an d (LN 1 ,N 2 Pk(LN 1 ,N 2 ))+ £ An u N 2 - (5-17) 

Proof Clearly P k (L NuN2 ) G ./W^. From f[5T2|) . ([531) and ([5TT3|> it follows 
that q k {x) G Aj\rx,jV2 f° r every k > N. This combined with equations ([5.15P 
and (|5.3ap shows that Pfc(L/Vi,Ar 2 )+ £ •^•JVi,Ar 2 - Similar argument gives the 
second statement in ([5. 17|) . Indeed, using that 

X + 2) G " 4 ' JVl,JV2 for every k ~ N ' 

and (|5.13p one can deduce that (x + l/x)q k (x) G Ajy^jVa f° r an k > N, which 
combined with uP k {u) = (x + l/x)q k (x) + (-l)*^ 1 ) + gives the 

second part of (|5.17j) . □ 

5.2. Characterization of LN lt N 2 m terms of the Toda flows. In this 
subsection we prove that the operators L^ lt N 2 can be characterized by the 
property (|5.17p . or equivalently by the vanishing of an appropriate linear 
combination of the Toda flows after a particular point. This is a discrete 
analog of the well-known fact that the rational solutions of the Korteweg-de 
Vries hierarchy are precisely the second-order differential operators, which 
are stationary under the KdV flows, after a particular point. 

Theorem 5.4. Let L = E + 6 n Id + a n E~ 1 be a nonconstant properly bor- 
dered second-order difference operator. Then the following conditions are 
equivalent. 

(i) The operator L can be obtained by a sequence of Darboux transfor- 
mations from the operator Lq = E + E~ l at the end points ±2 of the 
spectrum, i.e. L = Ljvi,at 2 f or some A 7 !, A^ G No with Aq + A2 > 
and for specific values of the free parameters in the Darboux process 

(EH). 

(ii) There exists N G N such that 

[P k (L) + ,L] = [(LP k (L))+,L] = for every k>N, (5.18) 
where P k (u) is the polynomial defined by (|5.16p . 



Proof. The implication (i)=>(ii) follows immediately from Proposition! 
Assume now that (15. 18f) holds for every k > N. Then L belongs to a rank-one 
commutative ring of difference operators Al and P k (L) + , (LP k (L)) + G Al 
for k > N. Let ^> n be the Baker function and X = Spec(.A) U {Qto^Qoo} 
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be the complete curve. Then, there exists a function f(P) on X with simple 
poles at and such that 

m> n {P) = f(P)*„(P). (5.19) 

Let us pick local parameters x~ l and y~ l near and Q^, respectively, 
such that 

( , p , _ J x + 1/x when P is in a neighborhood of , s 

| y + l/y when P is in a neighborhood of Qoo- 

Then for each k > N there exist functions f2fc+i(P) an d f2fc+2(-P) with poles 
only at such that 

P fc (L) + 1< n (P) = hk+i(P)*n(P) (5.21a) 
( J LP ifc (L)) + * n (P) = f 2fc+2 (P)*„(P). (5.21b) 

From (|5.20p and (|5. 15|) it follows that near and these functions have 
the following expansions: 

f2 f fe+ i(^) = %(^) + Ei ( 5 - 22a ) 



1 XJ 



OO c/fc 



i=o 



and 



f+ +2 (x) = (z + I) ft (x) + (-l) fc ( 2 * + ') + E I ( 5 - 23a ) 



oo rk 

f2k +2 (y) = E ( 5 - 23b ) 

j=0 y 

Let us consider now the function 

Bfc = fafe+if ~~ f2fc+2 G j4l- 

A straightforward computation using (|5.20p . (|5.22p and (|5.23p shows that 
near and we have 

fljf (*) = «f + ("I)" ( 2 \ + X ) + 0(l/x) (5.24a) 

0* GO = ^ + " T'o) + 0(l/y). (5.24b) 

If 7^ 0, the function will correspond to a difference operator from 
Al with support [—1,0], which would imply that L is a constant coefficient 
operator, contrary to our assumption. Thus 

5'q = (5.25) 
and Qk(P) is a constant (depending on k). Next, we define 

h = hk + s - (f - 4)f 2fc+ i + ("l) fe ( 2 \ + *)f e A L . 
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We have 



t } +(x) = -8 k x-8 2 * + 0(l/x) (5.26a) 

r k {y) = (s"l + (-i) fc ( 2 \ + X )) y - 4 + o(i/y), (5.26b) 



and the same argument as above shows that f)fc + 5 k f + b\ = 0, i.e. 



h, +:. ^ (f 2 - 4)f 2fc+ i + (-l) fc ( 1 ) f + 5 k f + 6l = 0. ( 5.27 > 



In particular, this equality implies that 



5N^ + (-l) fc+1 ( 2 \ +1 J (5-28) 
5 k = 6% (5.29) 



We want to show now that for every k > N, around we have 



oo 

w»)=«*( l /»)+E^ ( 5 - 3 °) 

From equations (|5.14p , (|5.20p , (|5.22ap and (|5.27p we deduce that coefficients 
5j in the expansion of the function f^. +1 (^) — %{x) = Y^jLi^j satisfy 
the following recurrence relations 

5 k = 5 k+1 + 35f (5.31a) 
S k +2 = S k+1 + 25 k - 8 k _ 2 for j > 2. (5.31b) 

The same argument shows that the same relations will be satisfied by the 
coefficients in the expansion of the function j 2k+l (y) — qk(l/y)- Thus, to 
prove that (|5,30p holds, it is enough to show that the coefficient of \/x 3 in 
Uk+i( x ) ~ 1k(x) is equal to the coefficient of l/y j in f 2k+1 (y) - QkO-fy) for 
j = 0, 1, 2 and every k > N. It is easy to see that the equality of these three 
coefficients is equivalent to equations (|5.25p . (|5.28p and (|5.29j) . completing 
the proof of (QUI) . 
Next, we denote 



p fc = 2f 2fc+ i - Pfc(f) 6 A L . 
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Using equations (|5.13j) . (|5.15|) . (j5.22aj> and (j5.30j) we see that pk(P) has the 
following expansions near and 

oo 

p+(x)=q k (x)-q k (l/x)+2j2 + 



/ -, \ 2k+l oo fik 



OO ^fc 
-. \ 2fe+l oo ^fc 

f--) +2 E^- ( 5 - 32b ) 

y) fri y 3 

We now use the fact that pk and f satisfy an algebraic relation of the form 

Pfc=PfcTl(f)+T 2 (f), (5.33) 

for some polynomials 71 and Tz- It is easy to show that 71(f) = 0. Indeed, 
if we assume that 71(f) = rf s + • • • for some nonzero constant r, then the 
function p&71(f) has the following expansions near and 



rx 2k+1+s + 0(x 2k+s ) in a neighborhood of Q 



J rx 2k+1+s + 0{x 2k+s ) 
PfcTx(f) = | _ ry 2k+i+s + 0{y 2k +s) 



+ 0(y ) in a neighborhood of Q r 

On the other hand, from equations (|5.32|) it is clear that if 

in a neighborhood of 
—ry + 0(y~ 1 ) in a neighborhood of Q^, 



Pk-T 2 (f) 



for some nonzero constant r, then I < 2k, leading to a contradiction. Thus 
71 (f ) = and therefore (|5.33p reduces to 

p|=T 2 (f). (5.34) 

We can show now that 5 k = for all j 6 N which implies that 72(f) = 
(f 2 — 4) 2fc+1 . Indeed, let S k Q ^ for some jo £ N and let jo be the minimal 
possible. Then we can rewrite (|5.34p as 

P 2 _ (f _ 4 )2fc+l = T2 ( f ) _ ( f 2 _ 4) 2fc+l_ (5 35) 

The left-hand side of (I5.35|) has the following expansions near and Q~ 

c 2_, f2 _ 4)2fc+ i = / 4^ x 2fc+1 "^ + 0(s 2fc -*>) in a neighborhood of Q+ 
Pfc 11 j \ -45 k y 2k+1 - 30 +0(y 2k ~ 10 ) in a neighborhood of Q~ . 

Again we get a contradiction because clearly the highest coefficients in the 
expansion of 7b (f) — (f 2 — 4) 2fc+1 near Q+ and Q~ must be equal. Thus 
(|5.34p becomes 

p 2 = (f 2_ 4)2fe+1 _ (5 _ 36) 
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Notice that the curve given by the last equation is exactly the spectral 
curve Spec(Ak,k) in I15.4|) . Thus, we can conclude that the spectral curve X 
is rational and choosing an appropriate parametrization we have 

A N;N C A L c C[x,x -1 ]. 

From Corollary 15.21 it follows that A = AN lt N 2 f° r some Nx,N2 G No with 
Ni + N2 > (because L is nonconstant). Moreover, as we saw at the 
beginning of the section, equation (|5.7p must hold. Let ^f n (x) be the Baker 
function for An u n 2 , i.e. the function in formulas (|5.3ap - (|5.3b|) . Then, the 
Baker function for Al is ^ n ( x ) = c~ n ^ n {x). Thus 

L#„(s) = f ~ (x + ^ + ^ * ft (a?). 

Equation (|5.36p shows that for k large enough the function 

. 2 I 2fc+l 

1 ( IN A , 

+ d -4 



must be the square of a rational function. It is easy to see that this can 
only happen when d = and c = ±1, which means that L = Ljv x jv a or 
L = Ln 2 ,N\ finishing the proof. □ 

Remark 5.5. We can easily reformulate Theorem [53] in terms of the vector 
fields Xfc of the Toda lattice hierarchy (|2.ip . For every k £ N let 

{0 if k is odd . n _. 

1 -ft- (5 - 37) 
1 it k is even, 

and let us define 

Yfc(L) -(L^-l)/2j)! ^ il (*-e fc -2j)! X ^' (L) ' (5 - 38) 

where Xj(L) = [(L J ) + ,L]. Then the operators Ln 1 ,n 2 with iVi + A^2 > 
can be characterized as the only nonconstant properly bordered second- 
order difference operators of the form L = E + 6 n Id + a n E~ x satisfying the 
constraints 

Yfe(L) = for every k large enough. (5.39) 

We end this section by giving an explicit formula for in terms of 
the vector fields X^- defined in Theorem I4.4L The proposition below will 
be needed in the next section when we want to characterize the operators 
having specific heat kernel expansions. 

Proposition 5.6. The vector fields Y^ defined by ()5.38p can be rewritten 
as a linear combination of the vector fields X' fc , corresponding to the flows of 
the system ()4.3|) . as follows 

Yfe = Yl (- 1 ) z ( 1 -f^J ( 5 - 4 °) 
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Proof. Using (H2J) and (IQHI) we get 

(k-e k )\ l(k ^l /2i ([(fc-l)/2]-j)! 
fe (L(fc-1)/2J)! ^ j! (fc-e fc -2j)! 

L(fc-i)/2j-i / k _ 2 j\ 

X S ( j ) X fc-2i-2j 



i=0 
L(fc-1)/2J 

1=0 

where 

(fc - £fc )! ^ (L(fc - l)/2j - j)!(fc - 2j)\ 

h ' 1 (L(fc-1)/2J)!^ i! 

We consider now separately the cases when k is odd or even. 

Case 1. Assume first that k is odd, i.e. k = 2s + 1. Then = and 
(|5.4ip gives 



(2^+1)1 (s-j) 



s! U i ! G-j)!(2s + l-/-j)! 
2a + l\ /-Z,- 2s -l + Z 

z J 2jF H -« ;1 

2 fl + l\(£+l-0, w/2a + l 

« / (-*)« 1 } \ i 

where in the last line we used the Chu-Vandermonde formula to evaluate 
the 2-Fi and (a); = a(o + 1) • • • (a + 1 — 1) denotes the shifted factorial. This 
completes the proof in the case when k is odd. 

Case 2. Assume now that k is even, i.e. k = 2s + 2 and therefore e k = 1. 
A similar computation as above gives 

(ga + l)! ' (-iy (s + l-j)! 
2S+2 '' *! J' ! a-i)!(2 S + 2-Z-i)! 

2S + 2 VF 1 ^- 2s - 2 + / ;l 



« 2ri \ -s-i 



2s + 2\(s + l- l)i _ , s - Z + 1 /2s + 2 



z ; (-s-i)j v y s + i 

which completes the proof. □ 

6. Finite heat kernel expansions 

In this section we prove that the operators Ln 1) n 2 constructed in the 
previous section can be characterized as the only operators for which the heat 
kernel can be written as a sum of only two Bessel functions with polynomial 
(in t) coefficients. 
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We work below with series of the form 

5>(r)/ r (2t), (6.1) 



r<k 



where p(r) is a polynomial in r, and before we state the main result, we 
establish several properties of the series f)6. 1 j) . Using (jl.7p we see that for 
r > we have 

tr ^ i*i 2) id 



' 3=0 J ' 



J*l 2 

Thus if |t| < T, then |J_ r (2i)| < ^ r -j-e' T ' 2 . Since for every polynomial 

p{r) and for every positive constant T the series YlrLo ~ l r ^ T r converges 
we conclude that the series (|6.ip converges absolutely and uniformly on 
bounded sets. 

The next lemma allows us to identify polynomials p(r) for which the sum 
of the series (|6.ip can be written in a simple closed form. 



Lemma 6.1. For every m 6 N and k £ Z we have 

4 m,fc 



t m I k (2t) = (-l) m Yl ^ih{2t), (6.2) 



i<k 

i=k+m (mod 2) 



where 



^= 4 * _ 1V n ^+i) 2 -^ 2 )- M 

|j|<m 
j=m (mod 2) 

Proof. The proof can be obtained by induction on m using the identity 

i(I fc _i(2t) - 4+i (2t)) = fcj fc (2t). (6.4) 

□ 

Notice that A™ ,k is an odd polynomial in i of degree 2m — 1. Thus, 
as an immediate application of the above lemma, we obtain the following 
corollary. 

Corollary 6.2. Let p(i) be a function which is an odd polynomial ini, when 
i is even or odd\j\ Then 

J2p(i)I i (2t)=p 1 (t)I k (2t)+p 2 (t)Ik-i(2t), 

i<k 

where pi(t) and P2(t) are polynomials of t such that pi(0) = P2(0) = 0. 

We can now formulate and prove the main result in this section. 

Theorem 6.3. Let L = E + 6 n Id + a n E~ l be a properly bordered second- 
order difference operator. Then the following conditions are equivalent. 



^Equivalently, we can say that p(i) = p(i) + ( — l) l p(i), where p(i) and p(i) are odd 
polynomials in i. 
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(i) The operator L can be obtained by a sequence of Darboux transfor- 
mations from the operator Lq = E + E^ 1 at the end points ±2 of 
the spectrum, i.e. L = LjVi,Ar 2 f or some Ni,N2 G No and for some 
specific values of the free parameters in the Darboux process (15. ip . 

(ii) The fundamental solution of the discrete heat equation (II. 2\ can be 
written as 

u(n,m;t) = (1 + pi(n, m; i))J n _ m (2i) +p 2 (™, m; i)I„._ m _i(2i), (6.5) 

where pi(n,m;t) andp2(n,m;t) are polynomials int with coefficients 
depending on n and m, such that pi(n,m;0) = p2(n,m;0) = 0. 

6.1. Finiteness of the heat kernel for the operators at 2 . The im- 
plication (i)=>(ii) is essentially proved in [15], except the fact that exactly 
two Bessel functions (I n -m and I ra _ m _i) are enough. This can be deduced 
from the arguments given there combined with Corollary 16.21 We explain 
below the main steps of the proof together with the essential ingredients 
from [EKE]. 

We consider the operator Ln u n 3 obtained by the sequence of Darboux 
transformations (|5.ip and the corresponding maximal commutative ring 
■An 1 ,n 2 °f difference operators that contains Ljvi,at 2 - From (|5.ip one can 
deduce that Ln x ,N2 an d Lq = E + E~ l are related by the intertwining rela- 
tion 

QLq = L^^Q, (6.6) 

where 

Q = Qn!+n 2 -i ■ ■ ■ QiQo- 

Equation (|6.6p implies that kevQ is preserved by Lq, i.e. Lo(ker Q) C kerQ. 
Conversely, one can show by induction that if two operators Lq and L = 
LN lt N 2 are related by (|6.6p for some difference operator Q, then L can be 
obtained by a sequence of Darboux transformations from Lq. The fact that 
we iterate the Darboux transformation only at ±2 means that the operator 
Lq restricted to ker Q has 2 and —2 as eigenvalues with multiplicities N± and 
N2, respectively. This allows us to reconstruct Q explicitly from its kernel 
as follows. We define functions <ftf, (j)^, ... , </>]$ and <f>J,<f>2,..., 4>Jj 2 such 
that 

{Lq - 2Id)0+ = for j = 1, 2, . . . , Nx 

(L + 2Id)0J = ij>J_ v for j = 1, 2, . . . , N 2 , 

with the convention that c/)q = 4>q = 0. Let Wr^ denote the discrete 
Wronskian (Casorati determinant) with respect to the variable n 

Wr A (/i, h, ...Jk) = det(A i - 1 / i ) 1 < 4J < fc . 
Then the operator Q, normalized to be monic, is defined by 
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The Baker function for the ring A-Ni n 2 can De written in terms of the 
operator Q by 

= ^i 7 4riF q(x " ) - <6 ' 8) 

Using the explicit form of the functions 4>^ , equations (|6.7p and (|6.8p one can 
show that ^ n (l/x) and x\]/* +1 (a;) are equal up to a multiplicative constant 
independent of x, i.e. 

* n (l/x) = C,x*; +1 (x). (6.9) 

Finally, we can use the above information to prove that \& n (x) and ^n( x ) 
satisfy the orthogonality relation 

^- £$> n (x)^* m+1 (x)dx = 5 ntm , (6.10) 

where C is a simple closed contour around the origin, avoiding the points 
x = ±1. The proof of (|6.10p can be obtained as follows. Notice first that 
^n(x) and x I / ^ +1 (x) have poles only at ±1. However the spectral curve has 
cusps at these points and since the differential * n (x)^'^ i+1 (x) dx is regular 
on the affine curve Spec(w4ATi,jv 2 ) we deduce that residues at x = ±1 are 
equal to zero. Using the explicit formulas for ^f n (x) and x I'^ +1 (x) we see 
that expanding around x = for m < n we have 

tf n (s)^ +1 (s) = - (<5„, m + O(x)) , 
x 

which establishes ()6.10|) for m < n. When m > n we replace x by 1/x in 
(I6.10P and applying (|6,9p we obtain zero by the bilinear identity (|2.9p . We 
refer the reader to [18] for detailed proofs of all statements in the above 
construction. 

From (|5.3ap and (|6.10p it follows that the fundamental solution for Ln 1 jy 2 
can be written as 

u(n,m;t) = -L / e*( x+ « i(x) dx. (6.11) 
2m J c 

Since 

e *(-+i) = ^x fc I fc (2t) 
we can deduce from (16. lip that 



u(n,m;t) = \ a k (n,m) 

± n—m—k 

where 



oo 
fc=0 



a fc (n,m) = ^ £ x fe - n+m * n (x)^ +1 (x) dx. (6.12) 

Using the last formula, we prove that for fc > 1, running over the even or 
the odd integers a^n, m) is an odd function of n — m — k with coefficients 
depending on n and m. The statement then will follow from Corollary 16.21 
The idea of the proof is to write x k ~ n+m in (I6.12|) as [ x k ~ n+m — p(x, n — 
m — k)] +p(x, n — m — k)iov appropriate p(x, n — m — k) which is a Laurent 
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polynomial of x, and an odd polynomial of n — m — k. We want to pick 
p(x,n — m—k) so that f(x) = x k ~ n+m — p{x, n—m — k) G An 1: n 2 - Then there 
exists a difference operator Lf = Ym=K- f Ji i( n )E l £ An%,N2 an d therefore 

k+ 

f(x)^ n (x) = L f ^ n (x)= ^ l^i(n)^n+i(x). 

l=K- 

If the interval [n + n + K + ] does not contain m we deduce from (|6,10p 
that 

^- £ f(x)V n (x)y* m+1 (x) dx = 0, 
and therefore formula (|6.12|) will give 

akin, m) = (p p(x, n — m — k)^ n (x)^* , i(x) dx, 

2m J c 

completing the proof. The main problem now is to construct a class of 
Laurent polynomials in A^^, which allows to implement the above idea. 
The key ingredient is the following proposition established in [15j . 

Proposition 6.4. Let N > max(iVi, N-z), and let Iq,1\,...,In be distinct 
nonzero integers, such that lj = Ik (mod 2) and lj + lk 7^ 0, for < j,k < N . 
Then 

N h 

S. — X 

Z> T~n 772 m G A N lt N 2 - (6-13) 

To complete the proof, we fix n and m, and we choose e = 1 or e = 2 so 
that k = e mod 2. Let N = max(iV"i, iVjj) and let us denote 

hS(j)= j ff j 2 -(j +21? 

JU) J + 2s JLX (j + 2s )2 - (J + 2/) 2 ' 

l^s 

Clearly, hj(j) is an odd polynomial of j. Now we can rewrite (|6.12|) as 
follows 

a k {n,m) = — I ( x k ~ n+m - V J#(ife - n + m)x J+2s \ ^ n {x)^>* x {x) dx 

+ £ ^; n + m) / - J+2 ^n(x)^ +1 (x) dx. (6.14) 
We need to define J depending only on n and m so that 

^- £ [x k ' n+m - p h s j(k - n + m)x J+2 ^ ^ n (x)^ +1 (x) dx = 0, 

(6-15) 

for every k, because then the right-hand side of equation (|6.14|) will clearly 
be an odd polynomial of n — m — k with coefficients depending on n and m. 
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Assume first that n < m and take J = m — n + e. If 1 < k < 2iV, then 
fc — n + m = J + 2s for some s G {0, 1, . . . , N — 1} and (|6.15j) is obvious 
because the polynomial 

N-i 

f(x) = x k-n+m _ h s j{k -n + m)x J+2s 

s=0 

is identically equal to 0. If k > 2N, then f(x) £ An 1: n 2 by Proposi- 
tion 16.41 and therefore f(x)^f n (x) can be written as a linear combination 
of {^k+m(x), ^k+m-l( x ): • • • > ^m+e(x)} and thus (|6.15p follows from (|6.10p . 

If n > m we can use a similar argument by taking J = n — m + e. This 
completes the proof of the implication (i)=>(ii). □ 

6.2. Characterization of Ljvi jy^ i n terms of the heat kernel. In this 
subsection we prove that (ii) implies (i) in Theorem 16.31 Assume first that 
L is a nonconstant second-order difference operator. The strategy of the 
proof is to show that (|5.39p holds using (|5,40|) and then apply Theorem 15.41 
and Remark 15.51 

From (|6.5p we see that 

u(n,n;t) = / (2t) + Pl (n, n; t)I (2t) +p 2 (n,n;t)I-. 1 (2t). (6.16) 

On the other hand, (|1.4p gives 

oo 

u{n, n; t) = I (2t) + ^ a k (n, n)J_ fe (2i). (6.17) 
fc=l 

Using equations (I6.16p - (l6.17p . the fact that the coefficients ak(n,n) in the 
expansion (16.170 are uniquely determined by u(n, n; t) and Lemma 16.11 we 
see that for k > 1, running over the even or the odd integers, ak{n,n) is 
an odd polynomial in k with coefficients depending on n. Similarly, using 
again (|6.5|) and writing u(n + 1, n; t) as 

u(n + 1, n; t) = I x {2t) + p 2 (n + 1, n; t)I (2t) + pi{n + 1, n; t)J_i(2t), 

we can conclude that for > 1, running over the even or the odd integers, 
a k+i{ n + Ijfi) is an odd polynomial in k with coefficients depending on n. 

Let 2N — 1 be the maximal degree of the four polynomials of k: n) 
when k is odd/even and ak+i(n + 1, n) when /c is odd/even. We show below 
that Y k (L) = for all k > 2N + 1. From (pTiOl) and (g^J) it follows that it 
suffices to prove that if f(x) is an odd polynomial of degree at most 2N — 1, 
then 

L(fc— 1)/2J 

E (-l) , (l-^ Jfc )(*)/(fc-20 = 0, (6.18) 

for every > 2 A r + 1 , where is defined by (I5.37P . 

If k is odd, i.e. k = 2s + 1, equation (|6.18p reduces to 

^(-l)^ 2s + 1 V( 2s + 1-2Z) = 0. (6.19) 
/=o ^ ' 
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Since / is odd, we see that the left-hand side of (|6.19l) is equal to 

2s+l 



\j2(-l) l ( 2S + 1 )f(2s + l-2l) 



and therefore (|6.19p is equivalent to 



E(-l) , ( 2 *^ 1 )/(2- + l-2/) = 0. 



(6.20) 



Equation (|6.20j) will follow immediately if we can show that for every poly- 
nomial F(x) of degree less that j, we have 



1=0 



E(-1)^)F(0 = 0. 



(6.21) 



The proof of the last identity is straightforward: 



D-d 



1=0 



l(3 
I 



F(l) 



1=0 



F{d z )Y,{-^{\y z 



2=0 



\z=0 



0. 



If k is even, i.e. k = 2s + 2, then (|6.18p is equivalent to 



E(-D 



1=0 



2s + 2 



I 



(2s + 2-2Z)/(2s + 2-2Z) 



2s+2 



z=o 



2s + 2 



(2s + 2-2Z)/(2s + 2-2/) = 0, 



which follows again from (|6.21|) . 

It remains to prove the statement when L has constant coefficients, i.e. 
L = E + 6Id + aE~ l , where a and b are constants and a ^ 0. In this case, 
we can write an explicit formula for u(n,m;t) and calculate m). For 

the fundamental solution we obtain 



u 



(n,m;t) = ^- <f> e(*+* +! :) t z nr ™- 1 dz, 



(6.22) 



c 



where C is a simple closed contour around the origin. Let w be such that 
z + b+ ^= w + -^ and w —* as z — > 0. Using (|6.22p . a short computation 
shows that 

' (1 — u> 2 )u>~ 1 ~ fc \ 



afc(0, 0) = res 



io=0 



(6.23) 



where 

q(w) = l-2bw + (2 + b 2 - 4a)w 2 - 2bw 3 + w 4 . 

In other words, a^(0,0) are the coefficients in the expansion of the function 
(1 — w 2 )l ' \J q(w) around w = 0. We know from Lemma 16.11 that if (|6.5|) 
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holds for m = n = then 0^(0,0) must be an odd polynomial in k for k 
odd or even. Let 



ajfe(0,0) 




when k is odd 
when k is even 



Then 



a k (0, 0)w k = fa(wd w ) w2j ~' + hiwd w ) ™ 2j 

k=l j=l j=l 

w w 2 polynomial of w 
= AH; o + P2{wd w )- 2 = n 2\K ' 

l — w z l — w z (l — w z y^ 

showing that YlkLi a k(®, 0)w k must be a rational function of w with de- 
nominator having zeros only at w = ±1. This implies that the only possible 
choices for q(w) in (16. 23ft are (1 — w) 4 , (1 + w) 4 and (1 — w 2 ) 2 . But the 
first two choices are clearly impossible (because «fc(0,0) become nonzero 
constants when k is even or odd), leading to c\(w) = (1 — w 2 ) 2 which is 
equivalent to 6 = 0, a = 1. Thus L = Lq completing the proof. □ 
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